EFFECTIVE MASS THEOREMS FOR NONLINEAR 
SCHRODINGER EQUATIONS 



CHRISTOF SPARBER 



Abstract. We consider time-dependent nonlinear Schrodinger equations sub- 
ject to smooth, lattice-periodic potentials plus additional confining potentials, 
slowly varying on the lattice scale. After an appropriate scaling we study the 
homogenization limit for vanishing lattice spacing. Assuming well prepared 
initial data, the resulting effective dynamics is governed by a homogenized 
nonlinear Schrodinger equation with an effective mass tensor depending on 
the initially chosen Bloch eigenvalue. The given results rigorously justify the 
use of the effective mass formalism for the description of Bose-Einstein con- 
densates on optical lattices. 
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1. Introduction 

Recent experiments on Bose-Einstein condensates (BECs) study the influence of 
optical lattices (or super-lattices) on the dynamics of the condensate, cf. ^1 El 
EHIES- The theoretical description of such systems is usually based on the famous 
Gross-Pitaevskii equation, i.e. 

j-2 

(1.1) ihdt4> = - — Ai) + V(xU + Uo(xU + Na M 2 'ilj, x 6 R 3 , iel, 
2m 

where m is the atomic mass, h is the Planck constant, N is the number of atoms 
in the condensate and 



(1.2) a 



with a £ R being the s-wave scattering length derived from the corresponding N- 
particle theory, cf. Depending on the sign of a, the condensate is said to 

be either repulsive (stable) or attractive (unstable). In the nonlinear Schrodinger 
equation (NLS) 1)1. l|l the potential Uq(x) models some given external confinement, 
whereas V(x) represents the lattice-potential, satisfying 

(1.3) V(x + j) = V(x), VxeR 3 , 7 er, 
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where T ~ Z 3 denotes some given regular lattice, generated through a basis {( 1 , C 2 j C3}) 
C, G M 3 , i.e. 

(1.4) £= j 7 GM 3 : 7 = ^7^, 7i ezj. 

Of course the nonlinear dynamics described by l|l.ll) can be highly involved. In the 
physics literature it is therefore frequently proposed, cf. |27l 1381 |4*T] , to consider 
the following simplifications: First it is assumed that the matter wave-field ip{t) 
can be characterized by a (fixed) central wave vector fco G M 3 and second one 
tries to capture the rapid oscillations in the wave function ^(t) by performing an 
asymptotic expansion in terms of Block waves Xn(y, kg) (see equation <|2.4|) below for 
their precise definition) . The center of mass of the wave function is then described 
by a slowly varying envelope function f(t,x), the dynamics of which is formally 
found to be governed by an effective-mass NLS. These types of approximations are 
well known in solid-state physics 8 , though mostly in a time-independent setting 
|35| . where one considers the motion of electrons in a crystal. It is the purpose 
of this work to rigorously justify the described approach in particular within the 
considered nonlinear context. 

To this end we shall consider a more general NLS than originally proposed in l|l.lf) . 
namely 

r h 2 

n _x I ihdti/j = - — Aip + V(x)ip + U (t,x)ip + a\ip\ 2a ?p, x G R d , iel, 

where a 6 R, and a G N. To motivate the choice a > 1, we note that for d < 
3 higher order nonlinearities are frequently used in the description of BECs, cf. 
[261 129| . Moreover different NLS type models of the form (|1.5fl also appear in 
nonlinear optics and laser physics, cf. 021 (see also [2] for a rigorous derivation). 
We assume ipj G L 2 (M. d ) to be normalized such that 



(1.6) / \ipi{x)\ z dx = 1. 

This normalization condition is henceforth conserved during the time-evolution. 
Again V is assumed to be periodic w.r.t. to some regular lattice F ~ Z d and 
Uq denotes some, in general time-dependent, smooth external potential. Now, we 
rescale the equation (|1.5fl . in order to precisely identify the asymptotic regime we 
shall be dealing with in the following. We have in mind a situation where the 
potential U$ is slowly varying on the lattice-scale corresponding to V. Hence, there 
are essentially two scales in this problem. First, the macroscopic length- and time- 
scale, denoted L and T respectively, which are introduced via Uq by rewriting it in 
the following dimensionlcss form 

(1.7) U (t,x) = ^U(t/T,x/L). 

In other words, the scaled potential U is such that a free particle of mass m under 
the influence of Uq will travel the distance L in the time unit T. On the other hand 
we can also introduce a couple of microscopic scales, A and r, via a rescaling of the 
periodic potential V such that 

(1.8) V{x) = ^V r (x/X). 
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The rescaled lattice T is henceforth generated through a basis {0}f=i) where Q 
Cj/\ and the microscopic time- unit r is then given by 

(1.9) 



h 

We consequently define two small dimensionless parameters, e and S respectively, 
as being the length- and time-ratios, i.e. 

(1-10) e=j, 5=1 

In the following, both of them are assumed to be small, i.e. s <C 1, S <C 1, but in 
general noi necessarily equal. Next we introduce new space- and time-variables x 
and t via 

(1-11) *=f, t=t 

and rescale the NLS 1|1.5|) in dimensionless form. Having in mind the normalization 
condition (|1.6(l we also need to rescale the wave function ip by 

(1.12) 4>{t,i) = L d l 2 il}{t,x). 

After multiplying (|1.5(l by T 2 /(mL 2 ), we consequently arrive at the following di- 
mensionless two-parameter model (where we again omit all " ~ " for simplicity): 

{ ihd t i> = - y + ^V r Q ^ + (t, a# + « 

Here we introduced two additional dimensionless parameters 

mL 2 mL da+2 
the former of which can be considered to be Planck's constant in the macroscopic 
variables. Note that the following important relation holds 

(1.15) e 2 = hS, 

connecting the ratio of the length-scales e with the corresponding time-scale ratio 
S. Finally, since we are aiming for nonlincarities of order 0(1) we shall impose from 
now on that 



lalT 2 mL da+2 

(1.16) \k\ = = 1, or, equivalently, T = J — — — , 

mL acr+z 1/ |a| 

hence relating the macroscopic length- and time-scales in a specific way. We remark 
that in the linear case a scaling analogous to (|1.13(l has been introduced in |37j . 

A brief discussion on several aspects of this scaling procedure is now in order: First 
note that if we choose h = e, hence, in view of (|1.15|) . e = 8, i.e. if we choose 
the same ratio for both, the length- and the time-scales, then the equation (|1.13|) 
simplifies to a one-parameter model given by 

2 

(1.17) isd t i)= ~Aip + Vr[^')p + U(t,x)ip + K\ip\ 2a '>p. 

This is nothing but the standard semi-classical scaling for (nonlinear) Schrodinger- 
type equations including an additional highly oscillatory periodic potential Vr- Re- 
cently the rigorous study of the corresponding asymptotic regime e — > 0, known as 
the combined semi- classical and adiabatic approximation, attracted lots of interest. 
In particular in the linear setting, i.e. n = 0, different mathematical approaches 
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are currently at hand, e.g., WKB-type expansions [HI 121) . Wigner transformation 
techniques |3()l I18j . and the so-called space- adiabatic perturbation theory |33l I44j . 
which gives the most sophisticated mathematical results so far. Including nonlin- 
ear effects, the literature is not so abundant. To the author's knowledge the only 
result in this direction is a recent paper by R. Carles, P. Markowich, and the author 
himself The results given there though are only valid for weak nonlinearities, 
in the sense that we need to assume k ~ O(e). Therefore a different rescaling of 
the original NLS l|1.5|l has been introduced in 

Remark 1.1. Additionally there exists a related work on the semi-classical limit of 
the so-called Schrodinger-Poisson system 5 in a crystal. There however additional 
assumptions have to be imposed which are out of the realm of the present setting 
(like truly mixed-state initial data). 

In the following our focus is not on the semi-classical regime though. Rather we 
shall study the asymptotic behavior of the scaled NLS (|1.13() for e <C 1 but with 
a fixed h of order 0(1). Note that, by (QBE , this implies S ~ 0(e 2 ), hence we 
are considering our system on a much larger macroscopic time-scale T than we 
do by fixing a macroscopic length-scale L via (| 1 . 1 0|> . In particular we are dealing 
with much larger times T as in the semi-classical studies described above. Roughly 
speaking the semi-classical regime can be seen to be an asymptotic description for 
ballistic scales, whereas we shall be dealing in the following with dispersive scales 
(sometimes also called diffusive scales). As we shall see, this indeed turns out to be 
the asymptotic regime where one can rigorously justify the effective-mass formalism 
discussed at the beginning. Remark however that, in contrast to what is noted in 
3 , the considered regime is not equivalent to the one obtained after rescaling time 
in the semi-classical equation (|1.17fl by t — > et. The reason being the different 
orders of magnitude in the external potential U and in the nonlinearity. 

To have a more concrete feeling of the involved time- and length-scales we come 
back to our original equation (|l.l|l . Thus we consider H1.13[l in d — 3, with a = 1 
and therefore k = Airh aT 2 / (mL 5 ) by (|1 . 14|> . A particular example for the periodic 
potentials used in physical experiments is then given by |14l I3ri| 

(1.18) V(x) = sin 2 (6x z ), &eR, 

' — ' m 

1=1 

where £ = (£i, £,2, £3) denotes the wave vector of the laser field which generates the 
optical lattice. Hence we readily identify A = (A 1; A2,A 3 ) as A; = 1/6- Moreover 
the slowly varying external potential Uq is usually modeled to be static and of 
harmonic oscillator type (isotropic or anisotropic), i.e. 

2 

(1.19) U a (x) = \x\ 2 , Lu eR,xeR 3 . 

In this natural choice for the macroscopic length -scale is therefore given 

by L = ao, where ao denotes the length of the harmonic oscillator ground state 
corresponding to Uq(x), i.e. 



(1.20) a 

V Loom 

The assumption s <C 1 is then of course equivalent to A C ao. In an actual physical 
experiment this requirement can be easily satisfied as a typical ground state length 
would be ao « 10~ 6 [m], whereas the wave vectors of the laser fields are usually 
tuned from 10 6 [l/m] to 10 9 [l/m], the latter case therefore being suitable in our 
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situation. The corresponding relation for the time-scales though is more subtle. 
From the condition (f 1 . 1 f>|> we see that T has to be chosen such 

T = * > t , since 6 < 1. 

47r7V|a|£ 

With t given by l|1.9fl this finally leads to the requirement 

»A 4 . 



47r7V|a| 

In particular, in the so-called moderate interaction regime, characterized by the fact 
that 4iTN\a\ « a @|, this is again equivalent to ao 3> A and in this case we compute 

a 1 



T 



AirN\a\u;l uj 2 ' 



Note that this is precisely what one would get in the corresponding linear situation. 

From a mathematical point of view the limit e — ► with S fixed, corresponds to the 
so-called homogenization limit of l|1.13fl . In view of the classical homogenization 
results as described in, e.g. |fill23|. the main new difficulty, apart from the appearing 
nonlinearity, stems from the large factor 1 je 2 in front of the periodic potential which 
furnishes a highly singularly perturbed term. It is therefore not a surprise that, 
even in the linear case, this type of homogenization problems have been rigorously 
studied only very recently |2]. In particular, (linear) time-dependent Schrodinger- 
type equations have been considered in j2j and in |37|. The latter result relies on the 
use of Wigner measures, a technique which can not be applied in the given nonlinear 
situation though. The former work is more closely related to ours, as it combines 
classical homogenization techniques, most notably two-scale convergence methods 
PP, with Bloch wave decomposition 13 . However, we want to stress the fact 
that the nonlinear case we shall be dealing with, is by no means a straightforward 
generalization of the linear results. More precisely, one should note that the scaling 
of (|1.13|) in general prohibits the derivation of suitable, i.e. uniformly in e, a-priori 
estimates, except for the basic L 2 estimate, corresponding to the conservation of 
mass. In the majority of cases, the derivation of such uniform estimates is crucial 
to gain sufficient control on the limiting behavior of the appearing nonlinearities, 
a problem which can not be handled by using weak-convergence methods (such as 
Wigner measures or two-scale convergence). 

Remark 1.2. Also, by the same reasons, our results do not fit in the framework 
of H-measures |43| . or G-convergence |34|. 

We note that in [21 EJ the authors propose the usage of a factorization principle 
in order to extend their results also to the nonlinear case. This approach though 
remains unproven there and moreover it is known to be applicable only in situations 
where the initial data ipi is concentrated at the minimum of the first Bloch band 
(see 0E] for more details on this). 

In comparison to that, the results given below are indeed independent of the number 
of the Bloch band and also they do not require ipi to be concentrated at a local 
minimum of the considered band. On the other hand we do need the initial data 
ipi to be well-prepared in a sense to be made more precise below, cf. Assumption 
14.31 Additionally we need to assume sufficient regularity on the potentials U, V 
as well as on the initial data ipj. The reason for these assumptions is on the 
one hand the fact that we shall use a more traditional multiple-scales expansion 
method, similar to those introduced in This approach will allow us to obtain, 



6 



C. SPARBER 



in a rather transparent way, an asymptotic description of ip(t) for small £ « 1, 
and also to determine the corresponding effective homogenized NLS. On the other 
hand, in order to prove that the given asymptotic solution is indeed stable under the 
nonlinear time-evolution governed by l|1.13fl . we shall adapt an approach originally 
introduced to prove the accuracy of nonlinear geometrical optics expansions, cf. 
|15ll2()ll3l?] for the most closely related results to the present work. The given proof 
will then again heavily rely on the fact that we have sufficient regularity properties 
and well prepared initial data. A similar strategy recently proved to be successful 
when applied to the weakly-nonlinear semi-classical situation studied in [Tl . The 
main goal of this paper though is not the introduction of new methods but rather 
a complete and correct treatment of the problem at hand. Moreover one should 
keep in mind that for completely arbitrary initial data ipi G L 2 (R d ), one can not 
expect an effective mass type dynamics to be valid. In other words, to obtain 
an equation of the form l|1.22fl the initial data ipi needs to be (asymptotically) of 
the same type as stated below, at least in leading order, and the additional well- 
preparedness assumptions we shall need, only concern the higher order terms within 
the asymptotic expansion. Remark that only in linear cases the superposition 
principle allows for more general states ipi. 

On the expense of not completely well defined assumptions (to be made precise 
later on) let us now state the typical effective mass theorem we shall prove in the 
following: 

Theorem 1.3. Let Vr and U be smooth, real-valued potentials, such that Vy is 
T -periodic and U is sub- quadratic. Assume that for some ko £ R, the initial data 
ipi is of the following form 

M*) = fi(*)Xn (J ,ko) e lk(rx/e + erf(x), 

where fj G iS(R d ;C), Xn(y,k) is an eigenfunction of Block's spectral cell problem, 
corresponding to a simple eigenvalue E n (k), and the corrector rf ~ 0(1) is such that 
ipi is well prepared, up to sufficiently large K G N, in the sense of Assumption \J~3[ 
below. Then there exists a eo > 0, such that for e < So the following asymptotic 
estimate holds 

(1.21) sup ||^)-«g(t)|| ia(B * ) =0(e) J r <r, 
te[-T ,T ] 

where t > is the maximum existence-time for a smooth solution f(t, x) for the 
homogenized NLS 

(1.22) ihd t f = ~ div(M*V)/ + U(x)4> + k* |/| 2CT /, 

with effective mass tensor M* = D 2 E n (k ) and an effective coupling k* G K, given 
by (|3.21(l below. Moreover, the leading order approximate solution Vq is found to be 

(1.23) v s (t,x) = f (t,x-fv k E n (k )j X n (j,*o) e *-/ £ e -^»( fc °)*A 2 . 

Remark that the estimate l|1.21[l implies a strong two-scale convergence statement 
as defined in Also note that, apart from the nonlinearity, our approach 

represents a refinement of the classical multiple-scales expansions given in , in the 
sense that we can include possibly appearing large drifts (clearly visible in (|1.23|l in 
the second argument of /) in order to resolve the underlying dispersive behavior. 
The possibility of large drifts present in the asymptotic solution can be seen as 
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the aftermath of the ballistic regime known from the semi-classical situation, a fact 
which has already been noticed in linear situations . 

The paper is now organized as follows: In Section[21we collect some preliminary re- 
sults and important notations used throughout this work. We then proceed similar 
to and first present in Section|3the multiple scales expansion method, whereas 
its nonlinear stability shall be proved in Section Q. 



2. Preliminaries 



(2.1) 



For simplicity we restrict ourselves in this work to static external potentials U = 
U(x), although all results could be generalized to the case of time-dependent po- 
tentials U(t,x) which are smooth w.r.t. t G R (indeed we could as well include 
smoothly time-dependent coupling factors n(t) G M). Thus we study in the follow- 
ing the asymptotic behavior as e — > of 

ihdti> = -yA^+ ^V r (J) V + U(x)if> + k H 2 >, 

All results given below are then valid for potentials which satisfy the following basic 
assumption. 

Assumption 2.1. The potentials U and Vr are such that i.e. Vr, U G C°°(R d ;R) ; 
and moreover they satisfy: 

(i) Vr is r -periodic: V T (x + 7) = V r (x), Vx G R d , 7 G T ~ Z d . 
(ii) U is sub-quadratic: d a U G L°°(R d ), \fa G N d , such that \a\>2. 

Remark 2.2. Remark that these are the same assumptions as used in JT]. In 
particular they include the cases of isotropic harmonic potentials U(x) = |x| 2 , as 
well as those corresponding to anisotropic ones, like U(x) — J2 UJ 'j x 'j- Moreover we 
can also take U to be identically zero, or include a linear component such as E ■ x, 
££K, modeling constant electric fields for example. 

We proceed by recalling Bloch's famous eigenvalue problem [7j. 

2.1. Bloch's eigenvalue problem. From now on we denote by € the elementary 
lattice cell, i.e. the centered fundamental domain of the lattice T, i.e. 
C d 



(2.2) C:={7£l": l = Y,7lCl, 



1=1 



1 1 

2' 2 



whereas the corresponding basic cell of the dual lattice will be denoted by (£*. In 
solid state physics €* is usually called the (first) Brillouin zone hence we shall also 
write 05 = €* . Also let us introduce the so-called Bloch Hamiltonian (or shifted 
Hamiltonian) given by 

1 

2 

Then Bloch's eigenvalue problem is given by the following spectral cell equation: 
' H r (k)xn(y, k) = E n {k) X n(y, k), n G N, y G £, 
Xn(y + l,k) = Xn(y, k), for 7 G r. 



(2.3) H r {k) := ^ (-iX7 y + k) 2 + V r (y), keR d . 



(2.4) 



and E n (k) G R, k G 03, is then called the n-th Bloch eigenvalue corresponding 
to the potential Vr- We shall now briefly collect some well known facts for this 
eigenvalue problem, cf. |31U40ll4*5] : 
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Since Vr is smooth and periodic, Hr{k), for every fixed k G 25, is self-adjoint on 
H 2 (T d ), T d = M. d /T, with compact resolvent. Hence its spectrum is given by 

a(H r (k)) = {E n {k) 6 K ; k e 25, n G N*}. 

The eigenvalues E n (k) can then be ordered according to their magnitude and mul- 
tiplicity, i.e. 

E x {k) < ... < E n (k) < E n+1 (k) < ... 

Moreover every E n (k) is periodic w.r.t. T* and it holds E n (k) = E n (—k). The set 

{E n (k) G K : E n (k) < E n+1 (k), fc G 25} 

is then usually named the nth-energy band (or Bloch band). The associated eigen- 
function, the so-called Bloch waves, Xn(v>k) form (for every fixed k G 25) a com- 
plete orthonormal basis in L 2 (€) and are smooth w.r.t. y G €. For the following 
we choose the usual normalization 



(2-5) (Xn(-,k),Xm(-,k)) L 2^ = J Xn(y,k)xm(y,k)dy = <S n , m , n, m G N. 

Regularity of the Xn w.r.t. their dependence on k G 25 is more subtle. It has 
been shown |3J that for any n G N there exists a closed subset 21 C 25 such that 
E n (k) is analytic. Similarly the eigenfunctions x„(-, k) are found to be analytic and 
quasi-periodic in k, for all k G D := 25\2l. Moreover it holds that 

(2.6) £„-i(fc) <£«(fc) <iWi (*), VfcGD. 

If this condition is satisfied for all k G 25 then E n (k) is said to be an isolated Bloch 
band |44j . Finally we remark that the set where one encounters the so-called band 
crossings, is indeed of measure zero, i.e. 

measSl = meas{fc G 25 | E n (k) = E m (k), n / m) = 0. 

Remark 2.3. Note that in the case d = 1 all band crossings can be removed 
through a proper analytic continuation of the bands, cf. |4()| . 

From the eigenvalue equation (|2.4|) we obtain the following useful identities: Dif- 
ferentiating the (|2.4() w.r.t. to k (assuming for the moment that everything is 
sufficiently smooth) yields 

(2.7) (V fc Hr(fc) - V k E n {k)) Xn + (H r (k) - E n (k))V kXn = 0. 

Hence, taking the scalar product with Xn, we obtain a formula for Vfe-E n (fc) by 

(Xn,^kH r (k)Xn) L 2 (e) = (Xn, {-i^y + k)Xn) L 2 {£) 
= V k E n (k), 

since H-p is self-adjoint. Similarly we obtain the following expression for the entries 
of the Hessian matrix D 2 E n (k) 

^ g ^ dl ]k E n {k) = 8jj + (xn, (-id yj + kj)xn) + (Xn, (~id m + h)d k] Xn) 
- (Xn, {d kl E n (k))d kj Xn + (d kj E n (k))d kl Xn), 

where Sjj denotes the Kronecker symbol for j, I — 1, . . . , d. (Below we shall assume 
E n (k) to be a simple eigenvalue which implies sufficient regularity to justify all 
differentiations above.) 
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2.2. Existence of smooth solutions for NLS. As a final preparatory step we 
state a basic existence and uniqueness result for NLS of the form (12.1(1 (see also 
[3 E2] for a general introduction) . 

Lemma 2.4. Let Assumvtion \2.1\ be satisfied, and let ipi G iS(M d ), the Schwartz 
space. Let s > d/2. Then there exists t — t(s,h) > and a unique solution 
V> G C(] - t,t[;H s (R d )) satisfying (|2~Tjl . Moreover, G C(] - i[; # s (K d )) /or 
any a G N d , seN, and the following conservation law holds: 

(2-10) | 11^)11^ =0. 

Proof. See the proof of Lemma 4.3 in □ 

Remark 2.5. In general, one cannot expect global- in-time existence for solutions 
to NLS. For example, if k is negative and if a > 2/d, finite time blow-up may occur, 
see e.g. ^302] ( see also 13 f° r the case a = 2/d). 



3. Multiple scales expansion 



We now establish the asymptotic behavior of ip(t), solution to 12. II for < £ < 1, 
by means of a multiple scales expansion. In the following h > is kept fixed, 
though we shall not simply put it equal to 1, since we want to keep track of its 
appearance in order to compare our results with the semi-classical situation of |llj , 
cf. Remark 13.91 below. Similar to we shall first consider the easier situation 
where no large drifts appear and then include them in a second step, using a more 
general asymptotic expansion method. 

3.1. Homogenization without drift. In this sub-section we seek an asymptotic 
expansion for solutions to i 1 . 1 3|) in the following form: 



(3.1) 



' / \ F ( x \ f f ko • x 
yj(t, x) = u yt, x, —J exp I i I h 

oo 

■J- e (t,x,y)~ 22e>Uj(t,x,y), 



where fco G M d is induced by the initial condition and (3 G M is some arbitrary 
constant to be determined below. The precise meaning of the symbol "~" in terms 
of an asymptotic series will be discussed in Section 0] below. Moreover we impose 
that u £ (t,x,y) G C satisfies 

u £ (v,y + 7 ) = u £ (;-,y), VyGM d , 7GT, 
in order to capture precisely those oscillations which are introduced via Vr- 

Remark 3.1. This particular form of a multiple-scale ansatz is suggested by the 
linear results given in [21 12) • Indeed one could have started with a more general 
ansatz, imposing appropriate periodicity or quasi-periodicity assumptions. It then 
turns out that one ends up with again the same form as given in (|3.1(l . Also, the 
ansatz (I3.1|l might not be so surprising when compared to the two-scale WKB- 
approach used in [rjl lllll?!] (see also the Remark 13 . 91 below) . 



As usual in asymptotic expansion methods we have to henceforth assume that the 
initial condition ipj is compatible with l|3.1|l . 
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Assumption 3.2 (Well-prepared initial data I). The initial data tpi is assumed to 
be of the following form: 

(3.2) ipi[x) = u e j (x, e lk °' x / £ , for some given k E R, 
where uf(x,y) is T-periodic w.r.t. to y and Uj E S(R 2d ) . 

Assuming for the moment that u e is sufficiently smooth we (formally) plug the 
ansatz (|3.1|l into i|2.1[l and compare equal powers in e. This yields 

(3.3) -^L u E + -L lU £ + Uu e + K\u e \ 2a u e = 0, 
e z e 

where the linear differential operators Lq and L\ are defined by 

Lqu £ := hf3 + h 2 H r (k ), 

(3.4) „ _ 

Liu 6 := -h, V x ■ V y u e + ih k ■ V x u e , 

with Hr(k) being the Bloch Hamiltonian as given in 1)2.3(1 . We also define 

h 2 

(3.5) L 2 u £ := -ih d t u e - — A x u £ + U(x)u e . 

Since u E ~ e 3 Uj we shall in the following expand equation l|3.3|l in powers of e 
and derive conditions on Uj, such that all resulting terms are zero up to sufficient 
high orders in e. 

Setting the leading order term, i.e. the term of order 0(e~ 2 ), equal to zero gives 

j3 

(3.6) H r {k )u + - "0 = 0, 

h 

from which we readily see that we need to choose 

(3.7) = -hE n (ko). 

Now, if we assume that E n (ko) is indeed a simple Bloch-eigenvalue then (|3.6I) 
implies that uq can be decomposed as 

(3.8) u (t,x,y)=f (t ) x)Xn(y,ko), V< el,i g R d , 

with some yet undetermined function f(t,x) E C. This now leads directly to the 
following important assumption 

Assumption 3.3 (Well-prepared initial data II). Initially, the leading order "am- 
plitude" uq is assumed to be concentrated in a single Bloch band E n (ko) correspond- 
ing to a simple eigenvalue of Hr(ko), i.e. 

(3.9) u (0, x, y) = / (0, x)\n(y, fco), 
where /o(0, •) = //(■) G 5(R d ;C) is some given initial data. 

An important consequence of E n (ko) being simple is that that in this case it is 
known to be infinitely diffcrentiable in a vicinity of feg. 

We have seen that in leading order ip(t) can be written as 

(3.10) X) ~ f (t, x) X n (J, fc ) e «o-»/« e -^(*o)*/e" + (e), 

where the /o is yet to be determined. To this end we proceed with our asymptotic 
expansion by setting terms of order 0{e~ l ) equal to zero. This yields 

(3.11) (H r {k ) - E n {ko))u\ = V x ■ V y u - ik ■ \7 x u , 



^That is, «| is smooth and rapidly decaying w.r.t. x and smooth w.r.t. y. 
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which by inserting ((3.8(1 can be rewritten as 

^ (H r {k Q ) - E n (k ))ui = - i^xfo ■ (-iVyXn + kaXn), 

= - iV^/o • Vfcif r (fco)Xn7 

where the second equality follows from the definition of Hp(k) ((2.3(1 . It remains to 
be asked if this equation is solvable for u\. By Fredholm's alternative the necessary 
and sufficient condition to do so is that the right hand side is orthogonal (in L 2 (£.)) 
to Xn{y, ko). Hence we have to impose that 

^ = - iV x fo ■ (Xn, VfeiJ r (fco)Xn} L 2 (C) 

= - i^xfo ■ V fe -E n (feo), 

where for the second equality we used the identity 1(2.8(1 . Thus we are lead to the 
restriction that kg has to be a critical point of E n (k), i.e. 

(3.14) V fc £„(fc ) = 0. 

This situation is analogous to the one discussed in the first part of (3J, though the 
arguments given there are different. Assuming that l|3.14|l indeed holds true, we get 
from 13.12fl . together with (|2.7|l . that the order 0(e)- corrector u\ in general can be 
written in the following form 

(3.15) ux(t,x,y) = -iV x f (t,x) ■ V k x n (y,k ) + fi(t,x)x n {y, k ), 

for any given function f\. Remark that we can not choose u\(0,x,y) completely 
arbitrary, once tto(0, x,y) is fixed, i.e. the initial data u e j given in Assumption 13.21 
(formally) has to be of the following form 

(3.16) uffoy) - (M0,x)+ef 1 (0,x))xn(y,k a )~ie\/ x f (0,x)-V kX n(y,k )+O(e 2 ), 

in order to be consistent with our asymptotic description. In the following we shall 
choose fi(0,x) = for simplicity. 

Proceeding with our e-expansion of l|3.3|l . we next consider terms of order 0(1) to 
obtain the following equation: 

(3.17) L u 2 + Liui + L 2 u + k\u \ 2,7 uq = 0. 
Again, by Fredholm's alternative, it can be solved for u 2 , iff 



(3.18) / Xn(y,k ) (Lxiix + L 2 u + K\uo\ 2cr uo) dy = 0. 

Plugging into this identity the precise forms of u$ and u\ , respectively defined in 
(13.8(1 and ((3.15(1 . and using formula 1(2.9(1 . given the fact that VfcS n (fco) = 0, we 
obtain after some lengthy but straightforward calculations the following solvability 
condition: 

(3 ig) | ihdtfo = ~y div a (M*V x )/ + U(x)f + k* |/ | 2ct / , 

I M t =o = -fr^)- 

This is nothing but the homogenized NLS, or the effective mass NLS, where the 
so-called effective mass tensor M* £ M. dKd is defined by 

(3.20) M* t :— d1 j kl E n {k ), j,l = l,...,d. 

Moreover the effective coupling constant k* 6 R within the n-th Bloch band is 
defined by 

(3.21) K*(k ) :=K f |Xn(y,fco)| 2CT+2 dy. 
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The effective NLS (|3.19|) describes the dispersive dynamics, as e — ► 0, of (|2.1|l 
for long macroscopic time-scales. However, it should not be confused with the so- 
called effective Hamiltonian as determined in |3XI I44| . Note that in general M* is 
neither positive nor definite, thus equation (|3.19|) in general also includes the class 
of so-called non-elliptic NLS [T71I231B2] . 

Remark 3.4. The formulas (|3.19() - (|3.21ll can be checked to be exactly the same as 
in the physics literature |381 l4Tj and moreover simplify to the ones given in [21 EH 
in the linear case. If M* is scalar its inverse m* — 1/M* is called the effective mass. 

In the next subsection we shall show how to get rid of the additional assumption 
(|3.14|l that kg is a critical point of E n (k). 



3.2. General situation including drifts. In order to generalize the expansion 
to situations where \7kE n (k) ^ we have to modify our multiple-scales ansatz. It 
turns out that instead of (|3.1(l we need to consider 

(3.22) ip{t,x) ~ u \t,x, -J exp U I — II, 

where u e ~ ^ s^Uj an( i the new spatial coordinate x is given by 

(3.23) x := x w(fco)t, with u)(ko) '■= VkE n (ko)- 

e 

Thus x comprises a macroscopically large drift with a drift-velocity proportional 
to Vfc£Vi(fco)- Note that the fast scale x/e remains unchanged, hence in situations 
where VkE n (k) = we are clearly back to our old situation. Similarly as before 
we plug H3.22|l into (|2.1(l which yields 

(3.24) -^L u e + -L lU e + L 2 u e + k\u e \ 2 ' j u e = 0. 
e z e 

where the linear differential operator Lo, L 2 arc defined as in (|3.4|l . I|3.5|) . respec- 
tively, but with x replaced by x. However instead of L\ we have 

(3.25) L x u e := -h 2 V £ ■ V y u e + ih 2 (k + uj(k )) ■ V^?/, 

Then, by exactly the same arguments as above we obtain that the leading order 
amplitude is given by 

(3.26) u (t,x,y)=fo(t ) x)xn(y,ko), Vfel,iel d . 
However, setting next the C(e _1 )-term equal to zero yields, instead of 13. 1211 . 

(3.27) (ffr(fco) - E n (k ))u! = -iV x f • (V fe i/ r (fc )Xr l - oj(k Q ) X n)- 
In this case, the solvability condition requires 

= (Xn, (V k H T (k ) - Uj(k ))Xn) L 2 {€) 

(3.28) 

= V fc £' n (fco) - w(fco), 

where we have used the normalization {\ n , Xn)L 2 (<i) = H3.28|) of course is 

identically fulfilled by definition of w(fco) : = V/ c £' Tl (fco). Thus, by identity (|2.7|) . we 
formally obtain the same C(e)-corrector ui 

(3.29) m(t,x,y) = -iV x f (t,x) ■ V k x n (y, h) + fi(t,x)x n {y, h), 

for any given function f\. (Like above we set /i to be identically zero at t = 0, for 
simplicity) Thus by transforming the x-coordinate into x by (|3.23(l we can now 
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proceed with our asymptotic expansion, having gained the additional freedom to 
include the case VkE n (ko) ^ 0. The equation of order 0(1) now gives 

(3.30) L u 2 + Lim + L 2 u + k\u \ 2,7 u = 0. 

It is clear now that as before, the corresponding solvability condition i.e. 

(3.31) J Xn(y, fco) (Ziui + L 2 u Q + k|u | 2o uo) dy = 0, 

yields a homogenized NLS equation of the same form as in l|3.19|l . namely 

(3.32) ihdtfo = - y div(M*V)/ + U(x)f + k* \f \ 2a fo- 

One can easily check that even though Vfe£ , „(/co) 7^ in this case, all additional 
terms appearing in Ij3.31|l . cancel out identically, hence equation (|3.32|) remains. 

Remark 3.5. In the physics literature |38l l4*T] the variable-transformation x — » 
x := x — huj{ko)t/E is sometimes reverted, leading to an additional convective 
term on the left hand side of (|3.32|) . This however can be considered only as a 
formal statement since consequently the large factor e -1 would appear then in the 
homogenized NLS, a somewhat inconsistent formalism. 

To prove the existence of smooth solutions for the homogenized NLS we shall impose 
the following ellipticity-assumption: 

Assumption 3.6 (Ellipticity). We assume that at fco £ 23 it holds: 

d 

(3.33) ^M*i = E d lkE n (k )^i > C\a\ 2 , for ^&R d ,C>0. 

k,l = l 

Clearly, condition l|3.33|l is valid if fco € © is indeed a local minimum of E n (k). It 
may very well be possible to relax the above assumption, cf. Remark 13.81 below. 
Here we mainly introduced it for definiteness, since under the condition H3.33JI it is 
then straightforward to prove that the effective NLS (|3~T§|) (or equivalently (f3~3l^ 
has a smooth solution, at least locally-in-time. 

Lemma 3.7. Let the Assurrwtions \2.l\ and \3.b\ be satisfied, and let fi 6 5(M d ). 
Then there exists t — r(h) > and a unique solution fo € C(] — r, t[; H s (M. d )), s > 
d/2, satisfying (|3~T3|) . or equivalently lj3~3l|l . Moreover, x a f £ C(] -r, r[; H s (R d )) 
for any a £ N d , s £ N. 

Proof. By Assumption Ij3.6|l we have that — div(M*V) is uniformly elliptic and 
since U is sub-quadratic the operator — div(M* V) + U is therefore well known to 
be essentially self-adjoint on C^ (R d ), cf. 40 . The existence of a smooth solution 
fo(t, •) £ H s (R d ), s > d/2, therefore follows by the same arguments as it does for 
the standard NLS [H]. The higher order regularity is then also proved similarly to, 
e.g., 01221 (see also the proof of Lemma 4.3 in for the main strategy). □ 

Remark 3.8. Indeed one can expect similar existence results to be true under much 
weaker conditions as given by Ij3.33|l . cf. 25 j. It is beyond the scope of this work 
though to study the weakest possible assumptions needed (including for example 
also degenerate cases) but rather we refer the reader to 021 and the references given 
therein. Here we only remark that in the case where M* is diagonal, the existence 
of smooth local-in-time solutions is known even in the non-elliptic case |17j . 
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Thus, at least for t G] — t, t[, we have that uo, and hence also ui, is smooth w.r.t. 
x, y and moreover in H s (M. d ) w.r.t. to x for any s <E N. It is clear that in general 
we can not expect smooth global-in-time solutions of the effective NLS l|3.19fl . A 
situation though were one indeed has globally smooth solutions, i.e. r = oo, is 
furnished by condition l|3.33|) together with assumption 12.11 and imposing that in 
addition k > 0. 

Remark 3.9. Let us briefly compare the obtained leading order asymptotic de- 
scription of ip(t) with the one derived in for the weakly nonlinear semi-classical 
scaling: Indeed if we formally set h = e in (|3.10() . 1|3.26[) . we obtain 

(3.34) t/,(t, x) ~ /„(*, x - u(k )t)xi (J, fco) e l ( fc o— ^(fco)*)/e + a ^ 

which is exactly of the same form as the two-scale WKB-ansatz used in JI] (see also 
[5ll21| for the linear case). In this case the highly oscillatory WKB-phase is simply 
given by <p(t,x) = fco • x — E n (ko)t. Note that this <$> solves the n-th band semi- 
classical Hamilton- J acobi equation with vanishing external field and plane wave 
initial data, i.e. 



(3.35) 



d t cj> + E n {V x 4>)=Q, 
fco ■ x. 



(3.36) 



k ^it=o 

Moreover with this choice of <f> (and since U vanishes) one easily checks that the 
transport equation for the leading order WKB-amplitude, as derived in |llj . sim- 
plifies to 

f d t fo + co(k ) ■ V x f Q = 0, 
\/o| t=0 = //(*). 
Clearly, the solution of (|3.36l) is then simply given by 

(3.37) fo(t,x) = f I (t,x-u(k )t), 
and hence consistent with our approach. 

3.3. Higher order expansions. We can henceforth proceed with our e-expansion 
of equation (|3.3|) . Denote the projector onto the nth-Bloch band corresponding to 
a simple eigenvalue E n (k) by 

Fn(fc) := \x n (y,k))(xn(y,k)\, 

(using the convenient Dirac notation) and moreover define 

Q n (fc) :=id-P„(fc). 

This operator is smooth in a vicinity of fco and hence, by elliptic inversion, a partial 
inverse for Lq = Lo(fco) can be defined on its range, i.e. L ~ 1 Q(fco) is well-defined, 
and smooth. Coming back to equation l|3.3Q(l we can decompose U2 as 

(3.38) u 2 (t,x,y) = f 2 (t,x)xn (y, fco)) + "2 (*i x > y), 
where the function f 2 is yet unknown and u 2 is such that 

P n (fco)^(t,a;,-) = (Xn(;k ),ui(t,x,-)) L 2 {£) = 0, V{t,x) e] -r,r[xM d . 
Now, U2 is determined by 13.30|l via 

(3.39) 11^2 = -L^Qniko) (L lUl + L 2 u Q + K\u \ 2a u )) , 

which implies G CQ — t, r[; H s (M. d )), since uq and u\ respectively are, by 
Lemma 13.71 As before, equation l|3.39|l henceforth induces a particular form of 
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the 0(e 2 )-corrector in the initial amplitude itf ~ ^2,e 3 Uj. The next higher or- 
der in e leads us to the following linear problem (after a Taylor-expansion of the 
nonlinearity around Mo): 

(3.40) L u 3 + L lU2 + L 2Ul + k ((2cr + 1)|u | 2ct Ui + 2cr|u | 2<T " 2 Uo«i) = 0. 

The corresponding solvability condition then determines f\(t,x) 6 C, i.e. the 
amplitude correpsonding to the polarized part of the first order amplitude u\(t, x, y) 
given in i|3.15[) . This then leads to a homogenized linear Schrodinger-type equation 
for fi(t, x), where we have the freedom to choose /i(0, x) = 0. 

By this procedure, all higher order terms Uj(t,x,y), j > 1, of the asymptotic 
solution (|3.22l) can be obtained and it is now clear that we can always choose 
gj(0,x), i.e. the polarized part of Uj(0,x,y), to be identically zero. In the globally 
periodic case, i.e. U(x) — and k = 0, the non-vanishing higher order terms 
Uj (t, x,y), j > 1, are found to be combinations of higher order derivatives w.r.t. k 
and x respectively, of Xn and /o, cf. [Hlll3j. Although this is no longer true in our 
case we still have that uj G C(] — r, r[; H s (R d )) for all j > 1. At each step though, 
an additional condition is imposed (recursively) for the initial data ipj . This can be 
seen analogous to the situation encountered in jl 1| and can be understood in the 
framework of so-called super- adiabatic subspaces as constructed in |33j . 

Remark 3.10. Note that the above given construction can be extended to the case 
where E n (k) is an m-fold degenerate family of eigenvalues, i.e. 



under the additional assumption that there exists a smooth orthonormal basis 
Xl{y^o)i ^ = l,...,m, of ranPj(fe), where Pj(fco) denotes the spectral projec- 
tor corresponding to E*(k). In this case the leading order asymptotic description 
would be 



(3.41) i>(t, x)~Y, /o,l(t, x) X i ( fc ) e *o-/e e -^„(feo)t/ £ 2 + 0(e)j 



In this case however we would be forced to consider a coupled system of homog- 
enized equations. The corresponding analysis is then analogous to the given one 
but requires rather tedious computations, a situation which we wanted to avoid for 
simplicity. For an extensive study of such situations in the linear case we refer to 
the last section of jHj- 



To prove that the above given multiple-scale expansion indeed yields a good approx- 
imation of the exact solution ip(t) for e <C 1, a nonlinear stability result is needed. 
Note that due to the scaling of l|2.1(l we can not hope for any uniform (w.r.t. e) 
bound in, say, H s (M. d ) for ip(t). On the other hand the uniform L 2 estimate (|2.1()|l 
is clearly not sufficient to pass to the limit in the nonlinearity. This motivates the 
introduction of the following e-scaled spaces: 

Definition 4.1. For s G N let 



E n (k ) = E4k ), Vne/cN, \I\ 



rn , 



n i 



4. Nonlinear stability of the asymptotic solution 




where 



II/ £ IIy/ : = E IM'W/ 6 



L 2 (B d ) ■ 



\a\ + \P\<3 
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Remark 4.2. Similar spaces, but without the extra weight e a have been used in 
the semi-classical study given in . Both variants can be seen as an extension of 
the i/^-spaces, defined by 



(4-1) ||/ E |U := E WWf 



L 2 (R d ) ■ 
\P\<a 

which in the context of geometrical optics expansion have been first introduced in 
PU| (see also and the references given therein). In our case the additional factor 
(ex) a is needed because we want to include sub-quadratic potentials U(x) (and due 
to our scaling we can not work in the A|-spaces introduced in If we would 

allow U(x) to grow only only sub- linearly we could as well work in H^. 

Notation. Let (a £ ) < e <i and (/3 e )o<e<i be two families of positive numbers. From 
now on we shall write 

* e < P E , 

if there exists a C > 0, independent of e s]0, 1] (but possibly dependent on other 
parameters), such that 

a £ <C(3 £ , for all e e] 0,1]. 

Since for the following results the value of h, appearing in (|2.1[1 . is indeed irrelevant 
we shall set h = 1 throughout this section. Moreover, we shall no longer distinguish 
between the usual spatial coordinate x and its shifted value x, since our results apply 
in both situations. 

Next, let us specify precisely the class of well-prepared initial data which we need 
to consider. 

Assumption 4.3 (Well-prepared initial data III). The initial data ipj satisfies 
Assumptions ^^ and \H.°A such that for some K £ N, it holds 

K 

(4.2) u e I {x) = Y i ^u j {x,y) " 1,Jl + l1 

3=0 

Moreover, with uq, u\ given by <|3.9fl . respectively, and with F(z) := |z| 2<T ; 

the functions Uj, j > 2, are recursively given by 



+ O (e + ) , in Y* for any seN. 

y=x/e 



7-2 , 3- 2 



-L 1 Q(k ) LxUj-i + L 2 Uj-2 + k - 2 F(u + s 



s=0 



After what we have encountered in the construction of higher order asymptotic 
solutions, this assumption should not come as a surprise. In the linear case the 
Assumption 14 . 31 is needed if one aims for refined asymptotic estimates. As we shall 
see, the inclusion of higher order asymptotics in our case is needed to control the 
nonlinear term in the proof of the stability result. To this end we need the following 
existence result for well prepared initial data: 

Lemma 4.4. There exists ipj G iS(M d ) such that Assumption \4.3\ holds true for any 
KeN. 

Proof. The proof follows from Borel's theorem, cf. Theorem 4.2 in [^j. □ 
For the following, define, the JV-th order asymptotic solution by 



(4.3) v £ N (t,x) := JV^ (t,x,^j e ik °- x 



/e e -ihE n (k )t/^ 



J=0 
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and moreover let 

(4.4) H e :=~A+^V r {^)+U(x) 

denote the linear part of the Hamiltonian operator. (Remark that the scaling of 
(14. 4|) is different from the standard semi-classical one as used in [111 I33| .) In the 
foregoing section we obtained the following preliminary result: 

Proposition 4.5. Let tpj satisfy Assumption \4- ^| for any K £ N and r > be the 

existence-time of smooth solutions to (|3.32[1 . Then for any N £ N, v N (t) solves 

(id t v N -H s v N = K ,\v N \ 2 °v N +e N r N , 

where r e N £ C(]-r,r[; H s (R d )), rj N+1 £ S(M d ) are such thatr% £ -r, r[; YJ>) 

and || r/ £ N+1 \\ ys = 0(1) /or any s £ N. 

The main result we shall prove is then given by the following theorem: 

Theorem 4.6. Let ipj satisfy Assumption \4-3\ for any K £ N, r > be the 
existence-time of smooth solutions to <|3.32[1 . and v £ N given by (|4.3[1 . Then for 
any To < t, there exists £o > such that for < e < Eq, the solution ip(t) to (|2.1|) 
is defined on the time-interval [— To, To] and moreover it holds 

(4.6) sup U(t) ~ v S N(t)\\ Y - =0(e N+1 ), 

tel-ToiTo] 

for any N £ N and seN. 

The above given results shows that if f(t) does not blow up in finite time, then 
neither does ip(t), at least for e sufficiently small. Further notice that if r = oo, then 
the above given estimate (|4.6|) holds for any bounded time- interval [t ,t ] £ M*, 
in contrast to the (nonlinear) semi-classical situation jl 1| where the appearance of 
caustics usually causes the WKB-approach to break down in finite time. Note that 
in this result, Assumption 14 . 31 on the initial data ipi has to be valid for any K £ N. 
We shall show in Proposition 14. 81 below how to relax this restriction. 



Proof. The proof is similar to the one given in Due to the different scaling of 
our equation, we shall present it here in more detail though, which moreover should 
benefit the reader. Define the difference between the exact and the asymptotic 
solution as 

w £ N (t,x) := ip(t,x) - v e N (t,x). 
Then, from (|2.1() and l|4.5() . wn solves 

/ id t w% = + k (IVf^ - \vn\ 2<t v n ) -e N r%, 

^ ' 1 e I N+l e 

{ w N\ t=0 = e Vn+i- 

By Lemma 13.71 and the well known Gagliardo-Nirenberg inequality, we have that 
v £ N is uniformly bounded in L°°([— To, To] x R d ). We shall prove now that is 
also bounded in L°°([—tq,to] x R d ), by using a continuity argument which shows 
that wnE is actually small in that space, for N sufficiently large. To this end we 
first note that the following important relation holds 

(4.8) \\n H > =e~ d/2 \\n Hi ze- d/2 \\r\\ Yr 
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where the scaling factor e d ^ 2 can be easily seen by Fourier transformation. This 
then directly leads us to an e-scaled Gagliardo-Nirenberg type inequality, i.e. 

(4-9) IMIl-cr-) £ lkllff-(R*) £ £- d/2 \W\v e * , ^ s > d/2, 

which we shall heavily use in the following. 

Multiplying the equation l|4.7|l by w e N , integrating over M. d , and taking the imaginary 
part yields 

(4.10) a t ||^(t)|| L2 < || - \v^v" N \\ L2 +e N \\r%(t)\\ L2 , 

since H e is self-adjoint and \k\ = 1 by (jl.ltjfl . To proceed further we recall the 
following Moser-type lemma, the proof of which is a straightforward generalization 
of those given in ^2 QUI : 

Lemma 4.7. Let R > 0. s E N, and F(z) = \z\ 2a z for a e N. Then there exists 
C = C(R, s, a, d) such that if v satisfies 

\\(ex) a (edfv\\ Lx(Rd) < R for all \a\ + \(3\ < s, 
and w satisfies \\w\\ Loa < R, then it holds 

\\(sxr(sdf(F(v + w)-F(v))\\ LHMd) <C WWWvWvw 

\a\ + \f3\<s \a\+\P\<s 

We shall now use this lemma to factor out w £ N in the right hand side of (|4.1()(l 
and then taking advantage of the smallness of the remainder. By construction, 
w%(0, x) — O (e Ar+1 ) in any Y* . By Lemma [3.71 we can find for fixed To < r an 
R > Q, such that if N + 1 > d/2, then 

(4.11) \\w%(t)\\ LOO < R, 

for e sufficiently small. Hence, as long as <|4.11ll holds, l|4.10(l and the above given 
Lemma 14.71 with s = 0, imply 

dt\\w%(t)\\ L2 <C\\w%(t)\\ L2 +Ce N \\r%(t)\\ L2 . 

Thus by a Gronwall type estimate, we get, as long as (|4.11|) holds, that 

(4-12) \\wm\\ L i<e N - 

for t < t. Next we shall show how to obtain similar estimates for the momenta 
and derivatives ofw £ N . Applying the operator eW x to l|4.7l) yields (where, as before, 
F{z) := \z\ 2<J z) 

id t (eV x w%) = H%e\7 x w e N ) + neV x {F{^) ~ F(v e N )) 
+ [eV x ,H £ ]w e N -e N+1 V x r%, 
and hence we obtain the following energy estimate 

d t \\eV x w%(t)\\ L2 <\\eV x (F^)~F(v%))\\ L2 + \\[sV x ,H £ }w%\\ L2 

On the other hand we compute from l|4.4|l that 

[eV x ,H e ] = ^xV T (j) +eV x U(x). 
Since VVr is bounded and Vt/ is sub-linear, 14.13fl consequently yields 

d t \\e\7 x w^ N (t)\\ L 2 < \\eV x (F(^) - F(v £ N ))\\ L2 + ^\\w%\\ L2 + \\sxw s N \\ L2 
+ e N \\eV x r%\\ L2 . 
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Thus, using again Lemma l4~7l (with s = 1) together with Proposition 14.51 and the 
estimate (|4.12|l we get 

(4.14) d t \\e\7 x w £ N (t)\\ L2 < ||eV x «^|| ia + \\exw e N \\ L2 + e N - 2 . 

(Remark that the difference in the last term as compared to the semi-classical 
estimate obtained in To obtain an estimate for j|exw| ir || i 2, we proceed anal- 

ogously to obtain the following moment estimate: 

d t \\exw%(t)\\ L2 < \\ex(F(iP)-F(v%))\\ L2 + \\lex,H £ }w%\\ L2 

+ e"\\exr £ N \\ L2 . 

But, since [ex,H £ ] = — eV x we get, as long as (|4.11|l holds 

d t \\exw £ N (t)\\ L2 < \\sx(F^)-F(v%))\\ L2 + \\eV x w £ N \\ L2 + e N \\exr%\\ L2 

<\\exw%(t)\\ L z + \\sV x w%\\ L2 +e N . 

Putting (|4.14|) and (|4.16[) together, we have 

dt (HeV^H^ + \\exw%(t)\\ L2 ) < \\eV x w s N \\ L * + \\exw%(t)\\ L2 + e N ~ 2 . 

Hence a Gronwall lemma yields 

(4-17) \\w £ N (t)\\ Y ^e N - 2 , 

as long as (|4.11|l holds and by induction one arrives at 

(4-18) \\w £ N (t)\\ Y s<e N - 2s . 

For s > d/2, and as long as (|4.11(l holds, the Gagliardo-Nirenberg type inequality 
(|4.9|) therefore implies 

Kr(*)ll*-CR<) ~ £ ~ d/2 W w N(t)\\ Y - Z e N - 2s - d / 2 . 

Hence if indeed N — 2s — d/2 > holds true, a continuity argument shows that 
(|4.11|l is valid up to times \t\ = t, provided e is sufficiently small. In particular, 
w e N , and hence ip, is well defined up to times \t\ — To < t, for < e < s(t). 

It remains to prove the estimate ()4.6JI . Fix s,N £ N and let si > s be such that 
si > d/2, as well as Ni > 2si + N + 1. From l|4.18jl we conclude that 

sup 11^)11^ <^- 2s ><£ N+1 . 

te[-T ,r ] 



Since iVi > N, it is therefore straightforward, that 

\Yf- 



SUp \\v E N (t)-V e Nl (t)\\ <r N +\ 



t€[— T ,To] 

and hence, we deduce that (|4.6(l holds for any s, N G N. □ 

In the above given proof, the initial data ipi is assumed to be well prepared up to 
any order K £ N. Indeed this rather strong assumption can be relaxed as the next 
result will show. To this end we introduce the following notation: 

Notation. For every a £ M we denote by [a] € N, the ceiling of a, i.e. the smallest 
integer which is larger than or equal to a. 

Proposition 4.8. Let ip(t) be the solution of l|2.1|l corresponding to an initial 
data tpi, which satisfies Assumption \4-S{ for any K G N. On the other hand, let 
ip{t) be the solution corresponding to an initial data ipj, where ipx * s such that 
Assumptions ^.^ is satisfied for K > [3d/2] . Then for any tq g] — r, t[, there exists 



20 



C. SPARBER 



£o > such that for < e < so, ip e {t) * s defined up to times \t\ < tq and moreover 
it holds: 



sup 

fe[— to,to] 



tp(t) - ip(t) y3 =0 (e K+1 - 2s ) , fors>0. 



Proof. Since the proof follows the lines of the one for Theorem 14.61 we shall be 
rather brief. Similarly as before we introduce 

w(t, x) := ip(t, x) — ip(t, x). 

Then w(t) solves 

id t w = H e w + k (\ijj\ 2a ip - \^\ 2a Pj , 

{y| t=Q = O (e K+1 ) in Y e s for anys e N. 

(Note that there is no remainder r £ N in this case.) We can then argue as in the 
above given proof. We have that initially it holds 

11-55(0, OIL- £ £~ d/2 |R0, Oily- < e K+1 - d/2 , provided s > |. 

With K + 1 > [3d/2] + 1 > d/2, the same arguments as in the proof of Theorem 14. 61 
yield 

as long as (|4.11|l holds. Since K + 1 > d, we can choose s > d/2 such that 
K + 1 - 2s > d/2, i.e. we can choose s such that K + 1 > [d/2 + 2s] = [3d/ 2}. 
Therefore the above given estimate and 1)4. 9|l show that 1 )4.110 holds up to times 
|i| = t 0) for e < 1. □ 

Theorem l4.6l and Proposition ^. 81 then finally lead to the following statement, prov- 
ing also Theorem ^21 

Corollary 4.9. If ipj satisfies Assumption \4.'^ with K > [3d/ 2], then there exists 
£o > such that for < e < £q, the solution ip(t) to 1)2.1)1 is defined on the 
time-interval [—To, To], for any tq < t and the following estimate holds: 

(4.19) sup ||^)-«o(*)II^(b*) =0(e). 

*£[— To,To] 

Additionally, if K > [3d/2], then we have 

(4.20) sup - vm\\L~(m*) =0{e). 
te[-T ,T ] 

Proof. From Proposition 14.81 we deduce that Theorem 14.61 holds with K > [3d/ 2] 
The L 2 estimate (|4.19|) then is nothing but ()4.6)l with N = s = 0. The L°° estimate 
(OU)l follows similarly from Theorem ED and (gTSJ. □ 

In other words we deduce that the solution of 1)2.1)1 can be, up to an error of 
0(e), approximated by the leading order asymptotic solution vq, obtained from 
the multiple scales expansion, if the initial data is well prepared, i.e. including 
correctors up to K — [3d/2], which is slightly stronger than what was required 
for the semi-classical result given in There the analogous condition for the 

correctors has been K > d. On the other hand one might guess that the leading 
order estimates 1)4.19)1 , (14.20)1 are true even if the initial data is only "correct" up to 
leading order. The techniques used in the above given proofs though, do not allow 
the conclusion that this is indeed the case. Note however, that the higher order 
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correctors in the initial data tend to zero, as e — ► 0, in practically every reasonable 
sense. 

Remark 4.10. Finally, let us remark that one could also study the semi-classical 
asymptotic behavior of the homogenized NLS, i.e. the limit h — » of l|3.19fl . 
although, in view of the given scaling arguments and Remark 13.91 the word "semi- 
classical" should rather be understood here in a purely mathematical sense. To 
this end, the well known WKB-type method derived in ^2] can be adapted under 
suitable conditions on M* . In this case however, one can only hope for local- in-time 
results, i.e. results up to caustics. The combined limit e/h — > 0, h — > though, 
seems to be more subtle, in particular due to the somewhat hidden dependence of 
r on h in the above given results. 

Acknowledgments. The author is grateful to R. Carles for helpful discussions on 
this work. 
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